g-thermostatistics and the analytical treatment 
of the ideal Fermi gas 

S. Martinez, F. Pennini, A. Plastino and M. Portesi 
Instituto de Fisica La Plata (IFLP), 
Universidad Nacional de La Plata (UNLP) and CONICET, 
C.C. 727, 1900 La Plata, Argentina 

Abstract 

We discuss relevant aspects of the exact g-thermostatistical treatment for an ideal Fermi sys- 
tem. The grand canonical exact generalized partition function is given for arbitrary values of 
the nonextensivity index q, and the ensuing statistics is derived. Special attention is paid to the 
mean occupation numbers of single-particle levels. Limiting instances of interest are discussed in 
some detail, namely, the thermodynamic limit, considering in particular both the high- and low- 
temperature regimes, and the approximate results pertaining to the case q ~ 1 (the conventional 
Fermi-Dirac statistics corresponds to q = 1). We compare our findings with previous Tsallis' 
literature. 
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I. INTRODUCTION 



Nonextensive thermostatistics [II 0, 0, 0, H E| constitutes a new paradigm for statistical 
mechanics. It is based on Tsallis' nonextensive information measure 



S q = k B l (1) 
g — 1 

where {p n } is a set of normalized probabilities and k B stands for Boltzmann constant (k B = 1 
hereafter). The real parameter q is called the index of nonextensivity, the conventional 
Boltzmann-Gibbs statistics being recovered in the limit q — > 1. 

The new theory comes in several flavors, though. Within the literature on Tsallis' ther- 
mostatistics, three possible choices are considered for the evaluation of expectation values in 
a nonextensive scenario. As a set of (nonextensive) expectation values are always regarded 
as constraints in the associated g-MaxEnt approach |8[, three different generalized proba- 
bility distributions ensue. Let p n (n = 1, . . . , W) stand for the microscopic probability that 
a system is in the n-th microstate, and consider the (classical) physical quantity O that in 
the microstate n adopts the value o n . The first choice 0] for the expectation value of O, 
used by Tsallis in his seminal paper, was the conventional one: Y^=iPn°n- The second 
choice 0, Y^=iPn °n, was regarded as the canonical definition until quite recently and is 
the only one that is guaranteed to yield, always, an analytical solution to the associated 
MaxEnt variational problem 0]; notice, however, that the average value of the identity op- 
erator is not equal to one. Elaborate studies of the so-called g-Fermi gas problem, which will 
constitute the focus of our attention here, have been performed using this "Curado-Tsallis 
flavor" [HQ El Qlli. Nowadays most authors consider that the third choice 0, , 



usually denoted as the Tsallis-Mendes-Plastino (TMP) one, 



sr^W q 

ln\ = ±^n=l Eli ° n (1} 
\V) q — w ~ \±) 

2-m'=l Pn' 

is the most appropriate definition. 

We employ the latter choice in order to accommodate the available a priori information 
and thus obtain the pertinent probability distribution via Jaynes' MaxEnt approach 

extremizing the g-entropy S q subject to normalization (^Y^ =l p n = lj and prior knowledge of 

a set of M nonextensive expectation values {(Oj) q , j = 1, . . . , M}. As usual the constrained 
extremization is accomplished by introducing M + 1 Lagrange multipliers; in practice two 
(equivalent) procedures can be followed to do this. First, the variational procedure followed 
by Tsallis-Mendes-Plastino in Ref. gives the Tsallis' probability distribution in the form 

fi/(i-?) 

Pn = ~ (3) 

where [0 

. . MEM^kdgjj ... 

In 1 w q \*) 

l^n'=\Pn' 

is the so-called configurational characteristic, and Z q = J2 n fn ^ represents a "pseudo" 
partition function (in the limit q — » 1 it goes to Z\ e^j= lXj ^ instead of Z\). f n should be 



2 



positive (otherwise f n = 0) in order to guarantee that probability p n be real for arbitrary 
q -Tsallis' cutoff condition [H, (l^-: as a consequence, the sum in Z q is restricted to those 
states for which /„ is positive. 

Notice that the expression obtained for p n following the TMP recipe is explicitly self- 
referential. This fact often leads to numerical difficulties in concrete applications (see, for 
instance, Ref. [US); more important, it obscures the underlying physics because the concomi- 
tant Lagrange multipliers loose their traditional physical meaning |2l| . Martinez et al. 2^ 



devised a way to circumvent these problems by recourse to the introduction of new, puta- 
tively optimal Lagrange multipliers (OLM) for the Tsallis' variational problem. The idea is 
to extremize the g-entropy with centered mean values (a legitimate alternative procedure) 
which entails recasting the constraints in the fashion 



w 



5>„ ff fo B -<Oi>,) = J = 1,-..,M (5) 

n=l 

The ensuing microscopic probabilities are formally given by Eq. Q, where now 

M 

/n = !-(!- {o ]n -(0 3 ) q ) (6) 
i=i 

In this way, the configurational characteristic obtained with the OLM recipe does not depend 
explicitly on the set of probabilities {p n }. 

It is obvious that the solution of a constrained extremizing problem via the celebrated 
Lagrange method depends exclusively on the functional form one is dealing with as well as 
on the constraints, the Lagrange multipliers being just auxiliary quantities to be eliminated 
at the end of the process. As a consequence, TMP and OLM results should coincide. 
However their manipulation is, in the latter instance, considerably simpler (notice that 
the OLM variational procedure (13, 0, ^2^, |25l |2(| |2tJ solves directly for the optimized 
Lagrange multipliers). Comparing then the TMP and OLM approaches for a given problem, 
one realizes that the resultant probabilities (as well as the pseudo partition functions) are 
identical if 

, (TMP) 

A, = ^y— q = Af MP) i = l,..,M (7) 



n=lP r ' 

where use has been made of the relation Y2 n Pn = Z q ~ q 0, I2| which is valid under the 
assumption of the knowledge available a priori. 



For the sake of completeness, let us write down [26[ the set of equations that constitute 
the basic information-theory relations in Jaynes' version of statistical mechanics 17, 



J-(lnZ,) = -<Cy, (9) 

for j = 1,...,M. Here the partition function Z q is defined by |22] lnZ q = \nZ q — 
^2jLi A? (Oj) q . While the above OLM equations involve ordinary logarithms, the analo- 



gous TMP relations |l6| employ the so-called g-logarithms, \n q x = (1 — x 1_9 )/(g — 1), thus 
in the latter case the basic Jaynes' relations are not recovered and the physical sense becomes 
somewhat obscured. 
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Our goal 



Motivated by the success of the OLM procedure, and in view of some recent and quite 
interesting applications of the quantum distributions (see, for instance, Ref. |28[ for anoma- 
lous behaviors in thermodynamic quantities for ideal Fermi gases below two dimensions), 
we wish to address here with such a technique the non-interacting Fermi-Dirac gas. The 
Tsallis generalized treatment of such system was originally advanced in Ref. using the 
Curado-Tsallis flavor. Biiyiikkihg et al. further investigated 12|| the generalized distribu- 
tion functions employing an approximation to deal with nonextensive quantum statistics 
called the Factorization Approach (FA). This approach, which comes out to give approxi- 
mate results in the region q ~ 1, is valid for a dilute gas ignoring the correlations between 
particles and regarding the states of different particles as statistically independent. Most 
succeeding works on the nonextensive treatment of quantum systems are based on these 
approximate generalized distribution functions. The nonextensive fermion distribution was 
further analyzed in Ref. in the context, again, of the second choice, but without recourse 
to that sort of approximations. It was seen there that the FA faces some difficulties. Other 
interesting studies on the subject have been presented by Ubriaco [l3[ and, quite recently, 
by Aragao-Rego et al. [2^ . employing the third-choice expectation values along the TMP 
lines. The latter work focuses attention upon the thermodynamic limit and provides elegant 
analytical results. 

The Fermi gas problem is rather cumbersome to treat using the TMP algorithm. We will 
show here that the OLM approach allows for an exact treatment of Fermi distributions in a 
nonextensive scenario, in a simpler way. This, in turn, will make it possible to re-discuss the 
nature of the approximation scheme of Biiyiikkihg et al. The paper is organized as follows: in 
Sec. |TT]we sketch the quantum version of the OLM technique. Sec. Illll is devoted specifically 
to the grand-canonical description of quantum gases in a nonextensive framework. Our main 
results concerning Fermi systems are developed in Sec. IIVI where a careful study of the mean 
occupancy of discrete single-particle levels is given. Next we present approximate results 
for values of q close to unity, and compare them with previous works on the subject. The 
thermodynamic limit is addressed in Sec. IV1| by appealing to integral transform methods 
that are summarized in the Appendix. Finally, some conclusions are drawn. 



II. THE OLM PROCEDURE IN QUANTUM LANGUAGE 

Since we are going to address the ideal Fermi gas, we have to adapt our nonextensive 
OLM-Tsallis statistical language to a quantal environment. Our main tool will be the equi- 
librium density operator p, that can be obtained by recourse to the Lagrange multipliers' 
method. Within the nonextensive framework one has to extremize the information mea- 
sure 0, S 

1 - Tr (pi) 

q-i 

subject to the normalization requirement and the assumed a priori knowledge of the gener- 
alized expectation values of, say M, relevant observables, namely 



S q [p\ = : r J ( 10 ) 



{6 ^ =T wf I 1 u 
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In the quantum version of the OLM instance the constraints are recast in the manner 

Tr(p) = 1 (12) 
= j = l,...,M (13) 



Tr 



p*(Oi-(Oi)< 



where the g-expectation values {(Oi) q , . . . , (Om) 9 } constitute the external a priori informa- 
tion. Performing the constrained extremization of Tsallis entropy one obtains 2^ 



1/(1-9) 



P 

where the quantal configurational characteristic has the form 

M 
3=1 

if the quantity in the right-hand side is positive definite, otherwise f q 



(14) 



(15) 



-cutoff condi- 



e quanti 

tion 0, 3]-. Here {Ai, . . . , Xm} stands for the set of optimal Lagrange multipliers, and 



we have defined for brevity the generalized deviation as 5 q O = O — (0} q . The normalizing 
factor in Eq. (|14|) corresponds to the OLM generalized partition function which is given, in 
analogy with the classical situation, by 



Tr 



f V(l-9) 
J a 



Tr 



M 

i=i 



(16) 



where the trace evaluation is to be performed with due caution in order to account for 
Tsallis' cutoff, and 

e q (x) = [1 + (1 - q)x} 1 ^ (17) 

is a generalization of the exponential function, which is recovered when q — > 1. Let us 
remark that the density operator a la OLM is not self-referential. 

It is to be pointed out that within the TMP framework one obtains from the normalization 
condition on the equilibrium density operator p the following relation that the OLM approach 
inherits B 0, 0, , namely, 



Tr 



f V(i-g) 

J a 



Tr 



f 9/(1-9) 

J a 



(18) 



Making use of this relation, one can obtain the value of the extremized g-entropy as 

S q = In, (Z q ) (19) 

For the sake of completeness, we can write down the generalized mean value of a quantum 
operator O in terms of the quantal configurational characteristic as 



Tr 


■^9/d-9) (J 


Tr 


f «/(!-«)" 

Jq 





(20) 



We recapitulate in the Appendix how to employ a quite useful method for calculating 
the generalized partition function and expectation values of relevant operators, by recourse 
to suitable integral representations. The procedure, which is based on the definition of the 
Euler gamma function, has been used in the literature by many authors 0, [i^, 27, 
3pl . l3lL I32L l33l | as it enables to express g-generalized quantities in terms of the conventional 
(q = 1) ones, thus providing an alternative analytic approach. 
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III. QUANTUM GAS IN A GENERALIZED GRAND CANONICAL ENSEMBLE 



It is our aim here that of developing formally the statistical description of quantum 
systems, particularly fermions, in a generalized framework. For this purpose, we will make 
use of the OLM-Tsallis version of nonextensive statistics. The Hamiltonian of the system is 
assumed to be of the form 

H = ^e k n k 

k 

while the number operator exhibits the appearance N = hk, where e k and denote, 
respectively, the energy and occupation number operator of the fc-th single-particle (s.p.) 
level for a discrete-energy spectrum, with e± < < ■ ■ ■ 

Following a grand-canonical-ensemble description, the configurational characteristic is 
given by 

f q = i - (1 - q)/3 (H - U q ) - (1 - q)a (N - N g ) (21) 

where the Lagrange multipliers (3 and a = —f3[i are related to the temperature and chemical 
potential /i of the system, respectively, and we have designated the generalized mean values 
of H and N by U q and N q , respectively. The OLM generalized grand partition function for 
this ideal quantum gas is obtained by inserting the last expression into Eq. 

In order to simplify the notation we now define e* k = e k — fi and the Legendre transform 
("free energy") H* = H — jj,N = ^2 k e* k hk- The quantities whose mean value is assumedly 
known, i.e. the internal energy and the particle number, are then combined in the fashion 

(H*) q = U q - fiN q = U* q (22) 

This notation allows one to treat the grand canonical ensemble as if it were the canonical one, 
but with grand canonical traces. In terms of the new ( "star" ) quantities, the configurational 
characteristic reads 

f q = i-(l-q)l3 (H*-U* q ) = ^g q (23) 

Pq 



where 



and 



Pq ~ i + (i - q)pu; (24) 



g q = i - (1 - q)(3 q H* (25) 



are auxiliary quantities. (It will be seen below that f3 q corresponds to the inverse 
temperature, while a q to the configurational characteristic, in a Curado-Tsallis treat- 



ment qShhEI-) 



Our partition function can therefore be written as 



a \ 1/(1-9) 

Z q =(^) TrfeVM) (26) 



= e q ({3U* q ) Tr [e q {-(5 q W 

where the last expression resembles the form of this function in conventional statistics. Once 
again, the trace entails the cutoff restriction which now reads: ((3 / j3 q ) g q should be positive 
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definite. From Eqs. (j2*U|l and (J23J), the scalar quantity [/* becomes 



u: 



Tr ( H* 



Tr 



?/(!-?) 



(27) 



Similar expressions hold for U q and N q separately, employing the corresponding operator 
inside the upper trace. A c omp arison of these generalized mean values with the concomitant 
Curado-Tsallis results 0, M, El allows the identification of the auxiliary quantity 
(3 q with the inverse temperature defined in the unnormalized context. 
The generalized heat capacity is defined as 



C 



dU n 



Vq 



dT 



d(3 



(28) 



After some algebra we can give it formally in the following fashion 



a 



Vq 



1 q(5 q [(g- 1 5 q H) q - (g- 1 5 q H6 q N} q (g^_ 5 q N) q / (g^_ {5 q Nf) q } 
l-ql- q p q [(g-i S q H) q - ($-1 6 q H 6 q N) q (g~^ S q N) q / (6 q N)*) q ] 



(29) 



IV. IDEAL FERMI GAS RESULTS 



In this section we present our fundamental results related to the nonextensive description 
of fermion systems, emphasizing the consequences of dealing with a discrete single-particle 
energy spectrum. We start by reminding that in the conventional statistical treatment of 



an ideal Fermi-Dirac (FD) gas, the grand partition function reads |3J, |3 



Z 1 = Tr 



-I3H* 



7V=0 {n k } 



(30) 



where the n^s take just two values (0 or 1) and the primed summation means that they 
add up to N. This fermionic partition function can be expanded as 



N=l fci<...<fcjv=l 



(31) 



in a form that, on the one hand, emphasizes the contribution of each possible value of N 
and, on the other, illustrates the manner of performing the trace operation. Note that, if 
the number of s.p. levels is finite, say K, both infinite sums in Eq. (j31|) nicely terminate 
when one reaches hi = K and N = K. This way of performing the pertinent sums, which 
is not the typical textbook procedure to deal with thermodynamic quantities in the case of 
Fermi systems, is of a general quantal {fermionic) character, no matter what the summands' 
content is. One could have, for instance 



Tr [p(-(3H*)] =l + J2 E V 

N=l fci<...<fcjv=l 



N 



(32) 
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involving an arbitrary analytical function <p. If ip is a generalized g-exponential we obtain 
Eq. ()33|) below. Indeed, in the framework of the g-thermostatistics one deals with the 
same sort of expansion, with a crucial difference: one faces, instead of the trace of ordinary 
exponentials, the trace of q- exponentials. Without recourse to the integral transform methods 
discussed in the Appendix, one can easily circumvent the main problem of the generalized 
Tsallis ' treatment: the fact that g-exponentials do not follow the distributive law with respect 
to sums over states, an important result of the present endeavor. 

Using the above considerations one thus recasts the Tsallis trace in Eq. (ffilj) in the fashion 



1/(1-9) 



E E 

N>1 ki<...<k N 



N 



1/(1-9)' 



(33) 



where, however, one has to take proper account of the cutoff requirement. Let us analyze 
this important point with some detail. 

The cutoff condition can be stated in the following manner 



q < 1 : 
q > 1 : 



N 

E 

i=i 

N 



i 



+ u* 



(34) 
(35) 



»=i 



for every possible (ordered) configuration [k\, . . . , fejy), and for all N = 1, . . . , K. In other 
words, those configurations not fulfilling the above inequalities do not contribute to the 
trace, i.e., the concomitant configurational characteristic is set to zero. Notice that the 
right side of both inequalities is nothing but [(1 — q)f3 q ]~ l and, in principle, could have 
positive or negative sign. Nevertheless, it can be easily seen that the requirement of positive 
probabilities -cutoff condition- is straightforwardly fulfilled by negative-definite "displaced" 
hamiltonians (whose spectrum is — fi) in the case q < 1 with U* > —[(1 — g)/3] _1 , and 
by positive-definite ones for g > 1 with U* < [(g — l)[3]~ l . In these two cases one can 
be sure that all configurations contribute with non-zero probability; but any other situation 
should be handled with some care. It is interesting to notice that the same sort of conditions 
are found in Ref. |33[ in the context of integral transform methods for the un-normalized 
(second flavor), canonical-ensemble problem. However, let us stress that in the present 
situation the analytic calculations can be fully implemented, with no other hardship than 
properly accounting for Tsallis cutoff as already mentioned, in summing over states in Z q 
(and other thermodynamic quantities). 

Let us now introduce, for the sake of brevity, the auxiliary scalar quantity 



N 



9 q {k 



1- 



k N ) = 1 - (1 - q)f3, 



i=l 



(36) 



which corresponds to the eigenvalue of the operator g q for the state with iV occupied levels: 
. . = rik N = 1, otherwise = 0. In terms of these new scalars the partition function 



reads 



PJ 



1/(1-9) 



ATM ki<...<k N 



(37) 
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Moreover, one can obtain U* from Eq. (|27|). evaluating the traces as exemplified with refer- 



ence to Z q . Thus 



TT* _ ^ ' " \ V (OO) 

i + E,> 1 E tl ,..<%#i--M' /N) 1 J 

Again, U q and iV 9 will also be given by similar expressions, replacing e£. by either e*^ or 
— fi, respectively. One easily ascertains that U* becomes = U\ — fiN in the extensive 
limit, and U* = — din Zi/d/3. Let us comment that in the case one wishes to determine, 
for instance, the dependence of the g-energy with temperature, one should consistently work 
out the above expressions in order to solve for the desired thermodynamic quantity. 

We tackle finally the evaluation of the g-mean value of the occupation number operator, 
(ni) q . To such an effect we, again, evaluate the traces in the form indicated above and obtain 
an exact expression for the generalized fermion occupation numbers 

_ Eiv>i J2 kl <...< kN (gjli faj g q (h, ...,k N ) «/(*-«) 
{ni)q ~ i + E*>i E kl< ... <kN 9 q (h,...,k N )^ (39) 

This expression can be further worked out employing a property of the auxiliary quantities 
g q (k\, . . . , /cat), which are not altered under a permutation of indices. Thus, we obtain 

, M _ 9,(1) ^ + E^> 2 E kl< ... <kN _ l{k ^ g q (kn ■ ■ ■ , j) ^ 
{Hl)q 1 + E kl< ... <kN 9g(ki,...,k N )"/^i) m 

Notice that (ni) q depends on the generalized internal energy U* through (3 q . 



V. APPROXIMATE RESULTS FOR q ~ 1: COMPARISON WITH PREVIOUS 
WORK 



Since the exact results discussed above exhibit a rather formidable appearance, and in 
order to gain a better grasp of the nonextensive thermostatistics of the Fermi gas, it is useful 
to consider the situation q — > 1. It is obligatory in this context to cite the pioneer work of 
Buyiikkihg et al. [121] . in which quantum gases are tackled in approximate fashion. They 
evaluate the partition function by recourse to the so-called Factorization Approach (FA), 
whose essential feature is that of ignoring interparticle correlations for the case of a dilute 
quantum gas. In other words, this is equivalent to treat the g-exponentials (|T7j) as if they 
were ordinary exponentials, what is approximately true when q ~ 1. The ensuing average 
occupation numbers have been widely employed in the literature It is to be stressed that 
these FA results were developed for the unnormalized second Tsallis-flavor mentioned in the 
introduction, namely, the Curado-Tsallis formulation (recently, some of us have brought 
this approximation up to date using the OLM recipes and applied the ensuing results to the 
black-body problem |27|). 

Since we are here in possession of exact results for the Fermi gas, we can indeed perform 
a check on the accuracy of the Factorization Approach. If we treat the g-exponentials as if 
they were ordinary exponentials we have 

JV 

g q (k u ...,k N ) VW.sJJ^) 

i=l 
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Strictly speaking, we are making use of an approximation for the g-exponential of a sum of 
entities in the form 

JV 

i=i 

which is valid for q sufficiently close to 1 such that the following inequality holds: 
|(1 — q) x if ~ Yl X T\ | •< 1. In our case, Xi = —j3 q e* k _. 

Therefore we obtain, from Eq. (|39|h the following FA- inspired approximate expression 




(fli) a « 77- ; (41) 

l + E Nil E h <...<,„^,1 1 [l-(l-«) ^ 3,4J s/,1 -" 

that can be cast as 



[i - a - rw 


{l + [l-(l-q)P g elf^} 








\ 



fa>* « — : . . . —1 ^ ( 42 ) 

leading straightforwardly to 



which is the Factorization Approach result of Ref. |12|, except that the power 1/(1 — q) in 
the result under an unnormalized context is changed to q/(l — q) under OLM strictures. It 
is clear then that the FA approximation is reasonably consistent in the q — > 1 limit. 

By using Eq. (|43J) one arrives to (formally) simple expressions for both the number of 
particles and the internal energy, 

N i « y 77T—^ ( 44 ) 

v [i - a - ?)/vr g/(1 ~ 9) + 1 

U w V - (45) 

9 r [i-(i-^^r ?/(H) +i 

similar to the ones obtained using the FA. We dare say that the present treatment is simpler 
than the one found in Ref. 111 ]. The simple appearance we are emphasizing here is deceptive, 
though. In order to perform any practical calculation one has to solve a coupled system due 
to the presence of (3 q . This problem, in turn, can be overcome by noticing that f3 q satisfies 
the approximate relation 

g-i+P-^E ^^V^ + i (46) 

which allows one to obtain (3 q in terms of f3 and, as a consequence, to decouple Eqs. (jUJ) 
and (gSJ). 

In this context, which we recall is valid for q ~ 1, we will now discuss the thermodynamic 
limit inspired in calculations performed by Ubriaco ^3|. We consider a system of massive 
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spinless particles in a volume V at temperature T. Going over to the thermodynamic limit 
in Eqs. (JUJ) and (@5J) we find that 

V ( 3 \ 3/2 

^^(^J (47) 

u A T ^Xff 2f ^ q) (48) 

where At = h/\/2wrnT is the usual thermal wavelength and z is the fugacity. We have 
introduced the following Fermi-like integral 



and 



x n ~ 1 



f*(z,o) = — — / s (49) 

JnV r(n)7 [i + ( ? _i) x _( g _i)|i n ^] 9 /( 9 -i) + i v ; 

Using the definition of Eq. (|17|h the first term in the denominator can be re-expressed as 
[e q (-x + (P q //3) lnz)]~ 9 so that 

i r°° x n ~ l 

f*( zq )= dx (50) 

Jn{ > q) T(n)J [e g (-x + (&//?) In z)p + l 1 ' 

and in this form it is easy to see that it gives the expected result z~ x e x when q — 1. 
The ratio (3/ (3 q can be obtained from the thermodynamic limit of Eq. (|4(-jj) as 



V 



3 



3\ 5/2 ( 3\ 3/2 



(51) 



Making use of the approximate results in Eqs. (J47|) and (J48j) we can express it in the fashion 

I „ l + (q-l)N q \nz 
q i + ( a _ i) N 2 

Then, we can write down the generalized energy per particle u q = U q /N q in a useful form as 

.,3 / 5 -/2fe9) l + (g-l)N,\nz 

'~~ 2 ^(^) 1 + (9 -i)iv,l|g ( ' 

that resembles the conventional result. 

Let us discuss the behavior of the generalized specific heat per particle in the present 
context. We have computed it following the usual procedure, by taking the temperature 
derivative of u q keeping the volume as well as N q fixed. The resulting expression -not given 
here- for Cv q /N q is somewhat involved; it can be written in terms of q — 1, N q , Inz and the 
Fermi-like integrals f*(z,q) with n = 1/2,3/2 and 5/2. An interesting study is to compare 
this generalized result against its conventional counterpart, Cy/N, to which it approaches 
when q — > 1. We have accomplished this comparison assuming that the value of the non- 
extensivity parameter q was close enough to 1 that we were allowed to use Taylor expansions 
for all the generalized quantities involved, up to first order in q — 1. We therefore defined 

C Vq /N q = C v /N [l + (q- 1) C (1) +0((q- l) 2 )] (54) 
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and obtained the relative first-order correction as a (rather complicated) function of z 
and N. In order to see the effects of non-extensivity on the specific heat for an ideal Fermi 
gas, we have considered the two extreme regimes of very low and very high temperatures. 
Our main conclusions are that: (i) Cy q (T = 0) = for arbitrary q; (ii) when T > (in which 
case Inz is of the order of /j,f/T ^> 1, where jxp stands for the Fermi energy), represents 
a positive contribution that behaves as iV Inz plus smaller terms, then Cv q /N q ^ Cy/N 
for q ^ 1; and (iii) when T — > oo (in which case z is approximately X^N/V <C 1), C^ 1 ' 
represents a negative contribution that also behaves mainly as iV Inz, then Cy q /N q ^ Cy/N 
for q 1 -in the case q > 1, the same finding is reported by Ubriaco jl3l] . 

Finally, let us point out the differences between the present results and the calculations 
of Ref. 13]. They manifest in the presence of powers of {(3/(3 q ) in Eqs. (|4"7|) . (|1SJ), and 
(I49J1 . These differences are ultimately due to the definition one chooses for the generalized 
mean values. As discussed in the previous sections, instead of using unnormalized mean 
values we work here within a normalized Tsallis framework, in its OLM version. Moreover, 
we compute, for consistency, the g-generalized expectation value for the number operator 
instead of dealing simply with (N)\. 



VI. THE THERMODYNAMIC LIMIT 

Let us discuss the thermodynamic limit in a more general context, and present our results 
for the internal energy and specific heat of an ideal fermion gas in the thermodynamic limit 
for any q > 1 -value (in the < q < 1 case we encounter a serious convergence problem on 
account of the cutoff condition. It is only easily tractable in the q — > 1 case). Following 
usual practice, in the limit of large volume (coordinate space) we can convert summations 
over discrete single-particle levels into integrations in phase space. If the energy spectrum 
of a particle in the gas is of the form e(p) = A |p| s with degeneracy g, we can write these 
integrals in terms continuous single-particle energies. In doing so we need the density of 
states, given by 

L d ? TT d / 2 1 / f \ d/s-i 
D(e)de = gt-^L—J- (±) de = ae b ~ l de (55) 
K ' y h d T(d/2) As \AJ V ; 

where we have assumed that the gas is contained in a hypercube of volume L d in a 
(i-dimensional space. In the case of massive spinless particles in 3D, one has a = 
2%V (2m) 3 ^ 2 / h 3 and b = 3/2; while for electrons in the ultra-relativistic limit one has 
a = 8irV/(hc) 3 and 6 = 3. The conventional (q = 1) grand partition function for an 
ideal fermion gas whose density of states is of the form , is given by [34J, |35j 

In Z 1 ({P, a}, V^oo) = aT(b) ^jf^- (56) 

with (3 > and b > 0. The function f n stands for the usual FD integral, which in terms of 
the fugacity z = e~ a reads 

i poo n— 1 °° yl 

fM =rh dx -i « + l = B" 1 )^ 
r (n) Jo z 1 e x + 1 ^ l n 

One can now write down the generalized grand partition function Z q making use of an 
integral representation (see Appendix). We consider here the real representation (Hilhorst 
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transform) of Eq. (jA3|) . To this end, the function Z\ is evaluated for the transformed 
Lagrange multipliers f3' = t(q — l)f3 and a' = t(q — l)a. Thus, in the thermodynamic 
limit we find, for any index q greater than one (and provided that the grand-canonical 
configurational characteristic is positive definite), the following exact expression 

1 r°° ( f ( p~ t ( t i~ l ) a ) \ 

2 ^7uaL "^"^(-t 1 -fa- w+ ° r w itt-iw ) < 57 > 

where I q stands for f3U q + aN q = PU*, from which one gets 1 — (q — l)I q = (3/ j3 q . 



A. Classical limit 

Let us consider now the case of sufficiently small values of z, which corresponds to the 
classical limit. In such a case, the FD integrals that appear in the different thermodynamic 
quantities of interest can be expanded in the fugacity, up to second order, as f n (z) ~ 
z — z 2 /2 n . The ensuing OLM generalized partition function becomes then 



Z q = e q (Q + a V e q (I q - a) ( \ _ ^ J 

a m r(^r-&) _ / I - (q - l)(Ig - 2a) \ 

(_L\ ^ 2a) { 2(q-l)(3 ) 



2 r 

+ a e q (I q -2a)^ ^—^ J +...(58) 

In the process we must take care of some restrictions on the parameters that arise from the 
use of the definition of the gamma function Eq. (jAljl . They are: i) q > 1, ii) 1 — (g— l)I q > 0, 
and iii) l/(q — 1) — 2b > 0. As a consequence, the present findings are valid within the region 
1 < q < min{l + 1/ (2b), 1 + 1/I q }- We can also obtain U q and in this case, in the form 
prescribed in Eq. (jA4|) . Keeping only terms up to first order in the fugacity, we are able to 
express the ratio between the generalized mean energy and particle number in the fashion 



U<1 ~N q - bl l-(g-l)6 lb9j 



from which we can obtain 



l + («-l)6(JV,-l) 1 

where one can easily recognize the appropriate result for the classical energy per particle 
when 5 = 1. In the present situation we also find I q = (b + a)N q /[l + (q — l)b(N q — 1)]. The 
implicit relation between (N) q and its corresponding Lagrange multiplier can be cast in the 
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following way 



(q - 1) a(N q - 1)\ / 1 + (q - 1) b(N q - 1 



i _ ( g _ i) ( aiVg + 6) ; V 1 - (? - !) - !) 

= aV(b) V- > e q (b) \^-J (61) 

from which one could, in principle, solve for a as a function of N q . Introducing this result 
into the expression given above for u q , one finally would obtain the generalized mean energy 
per particle as a function of temperature (and of N q ). Therefore, in this limit we are able 
to overcome the eventual possible complications introduced, within the OLM formalism, by 
the presence of the term I q = YljLi (Qj)q m ^ ne expressions for the generalized partition 
function and other thermodynamic quantities. 



B. Low temperature limit 

Let us face now the particular situation in which the fugacity is extremely large: 

z 1 i.e. a — > — oo 

which corresponds to the low temperature limit with finite chemical potential. Notice that, 
in this case, the transformed variable e~ t ^ q ~ 1 > a becomes very large as well. We can make use 
of Sommerfeld's lemma 134, 361 



fn(z) 



(In z) n L 7T 2 , . I „„ / 1 



r(n + i) 



1 + — n(n - 1)- + O 



6 (lnz) 2 \(lnz 



(62) 



which gives an asymptotic expansion of the FD integrals for z ^> 1. As can be seen, the 
dominant term in the expansion is of order (— a) n . We can evaluate the integral in Eq. (|57jl 
up to this order, getting an expression valid for all q > 1 in the low temperature limit: 



b 

-a 



z *- e *\h- {- ) (-) ] (63) 



Additional quantities evaluated up to the same order are 



U q ~— — fi b+1 and N q ~% b (64) 
o + l o 

independently of the value of q, in agreement with previous results obtained using the 
OLM formalism ji^. Using the above relations we can derive the Fermi potential as fip — 
(N g b/a) 1 / b ; and we can also identify the Fermi temperature as Tp = [if = (N q b/ a) l l b . 
Besides, it is easy to obtain I'f 1 and then simplify the expression for Z q: leading to 

Z q ^e q (2I^)=e q (- l( ^ T) (3^ (65) 
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It is worth pointing out that the mean energy per particle in the OLM-Tsallis framework 
does not depend explicitly on the value of q. Indeed, in the limit under consideration we 
have U^/Nq = n F b/(b + 1). 

Performing calculations up to the next order of approximation we arrive, after a little 
algebra, at the following results which are valid for any q > 1: 



TV ~ - u 
b 



71 



\ + — {b-l)b 
6 



T- 



+ . . . 



(66) 



where the quantity X = 1 - (q- l)/3{U q -fi[N q -afi b /b(b+ 1)]} = \-{q-l)(I q - 4° } ) is very 
close to unity. Actually, the second term in X vanishes if one keeps just terms corresponding 
to the lowest-order approximation, which eliminates the term I q introduced by the OLM 
procedure (see Eqs. ()13|h (JT3J), and (|7T|) ). For the chemical potential we find 



(2) 



71 

~6 



(6-1) 



T 



l-(ff-l) 



(67) 



whereas the generalized mean energy becomes 

= 



1 



6+1 



1 + 



I) 



T 



(68) 



From this expression we derive the generalized specific heat at constant volume 

,2 rp ] _ 



C v ^ = N a ^b^ 



3 T F 1 - {q - 1) 



(69) 



where we have assumed X 2 = 1. Notice the linear behavior of the specific heat with temper- 
ature and that the only difference with the conventional results (see, e.g., Ref. 34]]) comes 
through a factor 1/(2 — q) which goes to unity in the limit q — > 1 + . 



VII. CONCLUSIONS 

In this communication we have presented an exact statistical treatment for the ideal Fermi 
system in the generalized, nonextensive thermostatistics framework of the third Tsallis- 
flavor, the TMP one. Our main innovation is that of employing the OLM approach to 
nonextensivity (13, 111 HI HL IE 2?|, which allows one to obtain analytical results unavail- 



able if one uses other algorithms that revolve around the concept of Tsallis' entropy. 

Additionally, we have i) solved in exact fashion the Fermi-TMP equations, ii) introduced a 
method for evaluating traces that bypasses the use of the Gamma representation, iii) devised 
a rather simple treatment of q ~ 1 instances, and iv) studied interesting limiting situations. 
More specifically, the exact generalized partition function in the grand canonical ensemble 
has been given, and we derived the ensuing statistics for arbitrary positive values of the 
nonextensivity index q. Several limit instances of interest were here discussed in some 
detail: 

1. the thermodynamic limit, 
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2. the case q ~ 1 (g = 1 corresponds to the conventional Fermi-Dirac statistics), 

3. the low temperature regime, where we obtained results that are independent of the 
specific g-values. 

In writing down the generalized expectation value for the occupation number operator we 
were able to explicitly display the correlation among the occupations of different levels, which 
is a typical nonextensive effect. Indeed, the distinct mean occupation numbers are seen to 
disentangle from each other as one approaches the conventional q = 1 statistics. Finally, we 
discussed the limits of validity of the Factorization Approach of Biiyiikkihg et al. [l2T |. 
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APPENDIX A: INTEGRAL REPRESENTATIONS 

In this appendix we summarize a practical method for calculating the generalized parti- 
tion function (flljj) and expectation values (|2U|l of relevant observables, by recourse to integral 
representations based on the definition of the Euler gamma function. We start with the iden- 
tity 

/>oo 

/ dtt u - l e- tTI = r]- u T(iy) (Al) 
Jo 

for ^(77) > and > (see for instance Ref. j3?J, page 342). One can then write 

1 /*oo 

fW-ti = —/-^ / dt t^r- 1 e-*?* (A2) 

r (ft) Jo 

with the restrictions that i) f q be positive (which is indeed always complied with because of 
Tsallis' cutoff condition |l9|), and ii) q > 1. The usefulness of the transformation becomes 
evident, as the power- law form is converted into an exponential factor. 

Evaluating the trace of the above expression -see final comments below- and introducing 
the OLM quantal configurational characteristic, one arrives at an integral representation 
(Hilhorst transform [30j] ) for the OLM generalized partition function (Jlrjj) of index q > 1, in 
the form 

1 /*°° 

Z q = ' dtt^e-^-^^^Z^iK}) (A3) 

r (ft) Jo 

Here the integrand contains the corresponding conventional partition function Z\ evaluated 
for the set of transformed Lagrange multipliers {X'^t) = t(q — l)Aj, 1 = 1,... , M}. Notice 
that this transformation preserves the sign of the Lagrange parameters (this fact is used in 
the text). The generalized expectation values (fTTj) can be expressed, for q > 1, in terms of 
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integrals involving Z\ and its derivative with respect to the associated Lagrange multiplier. 
Indeed, one realizes that 



rdt t^ e -*[i+(i-?)E fe A fc (6 fc ) g ] dZl ({\i})/d\' 
/q.\ = _ £0 v L l} " J (A4) 

Let us mention that one can obtain alternative integral representations for Z q and {Oj) q 
in the range < q < 1 by recourse, for instance, to the following complex representation of 
the Euler gamma function 

dt (C + it)'" e {t:+it)v = 2tt n v - 1 /T(v) (A5) 
for rj > 0, > and jft(C ) > (|37|. p. 343). In this case (that may be called the 



Prato-Lenzi transform 31, 32]) the partition function can be written as 

Z q = —v( r dt (1 + ^)^i- 1 e( 1 + i ')[ 1 +( 1 -^S fc A fc (6 fc >, ] Zl ({Xi}) (A6) 
2vr \q-lj J-oo 

where Aj(t) = (1 + it)(l — q)\ for each i — 1, . . . , M, and the region of validity is q < 1 or 
q > 2. The associated mean values take the form 

j^dt (1 + it) e (i+«)[i+(i-9)E*A fc ((3 fc > 9 ] ^({Ai})/^- 
^ 9 ~ J^dt (1 +^)^ e ( 1 + 4 *)[ 1 +( 1 -^E fc A fc (o fc ),] ^({AJ) (A?) 

for j = 1, . . . , M and q < 1. 

Some additional remarks are necessary. A quite detailed analysis concer ning integral 
representations for the g-thermostatistics can be found in the recent preprint [23j by Solis 
and Esguerra, who pay special attention to practical details of the representations discussed 
above. One has to make sure that all states are contributing to the evaluation of the trace 
-in the sense that there is no cutoff- in order to get Z\ on the r.h.s. of Eqs. (jA3|) or (|A6|) . 
Solis and Esguerra point out that this fact has not been taken into account in the majority 
of Tsallis-related works (see, for instance, |29]). In our particular case, the condition to be 
imposed in order to ensure that the Hilhorst-type representation can be safely used can be 
stated in the following terms: f q as given by Eq. (}2*Tj) should be positive definite for all 
states, i.e., the lowest energy eigenvalue (or the greatest lower bound of the Hamiltonian) 
should be greater than or equal to ^(A^ — N q ) + U q — l/[fl(q— 1)]. (Notice that the simpler 
requirement given in Ref. [33j | . namely that H be greater than — 1/ [/?(<? — 1)], originates 
in the fact that a canonical-ensemble description is performed and also that unnormalized, 
Curado-Tsallis mean values are employed.) 
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